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ABSTRACT 

We have investigated the thermally induced proton/deuteron exchange in mixed amor¬ 
phous H 2 0:D 2 0 ices by monitoring the change in intensity of characteristic vibrational 
bending modes of H 2 0, HDO, and D 2 0 with time and as function of temperature. 
The experiments have been performed using an ultra-high vacuum setup equipped 
with an infrared spectrometer that is used to investigate the spectral evolution of ho¬ 
mogeneously mixed ice upon co-deposition in thin films, for temperatures in the 90 to 
140 K domain. With this non-energetic detection method we find a significantly lower 
activation energy for H/D exchange - 3840 ± 125 K - than previously reported. Very 
likely this is due to the amorphous nature of the interstellar ice analogues involved. 
This provides reactive timescales (r < 10 4 years at T > 70 K) fast enough for the 
process to be important in interstellar environments. Consequently, an astronomical 
detection of D 2 0 will be even more challenging because of its potential to react with 
H 2 0 to form HDO. Furthermore, additional experiments, along with previous studies, 
show that proton/deuteron swapping also occurs in ice mixtures of water with other 
hydrogen bonded molecules, in particular on the OH and NH moieties. We conclude 
that H/D exchange in ices is a more general process that should be incorporated into 
ice models that are applied to protoplanetary disks or to simulate the warming up of 
cometary ices in their passage of the perihelion, to examine the extent of its influence 
on the final deuteron over hydrogen ratio. 

Key words: astrochemistry - methods: laboratory: molecular - methods: laboratory: 
solid state - solid state: volatile - ISM: molecules 


1 INTRODUCTION 


The delivery of water to Earth is far from understood. Sev¬ 
eral hypotheses have been put forward, among which deliv¬ 
ery by comets and/or asteroids during the so-called ‘late ve¬ 
neer’ stage (Morbidelli et al. 2000 Cleeves et al. 2014), direct 


adsorption of water onto grains prior to planetary accretion 
(Muralidh aran et al.|2008 l, and in the recent ‘Grand Tack’ 
model, water delivery during the formation phase of terres¬ 
trial planets (O’Brien et al. |2014 1. Since the ratio between 
deuterons and protons incorporated in water molecules is 
elevated in the Earth Ocean with respect to that of the 
initial bulk solar composition, the D/H fraction in various 
molecules is a logical tracer for the proposed origins. This 
concerns particularly the ratio of HDO over H 2 0. The D 2 0 
abundance is expected to be low, because the amount of deu- 


* corresponding author 


terium present in the interstellar medium is a factor 100 000 
less than that of hydrogen. In fact the number of astronomi¬ 
cal D 2 0 detections is very limited: in the cold envelope layer 
surrounding solar-type protostar IRAS 16293-2422 (Butner 


et al.||2007 Vastel et al. 2010 Coutens et al.|[2013 1 and in 

the warmer regions surrounding the Class 0 protostar NGC 
1333 IRAS2A ( |Coutens et al.|2014 i. 


It should be noted, though, that a direct comparison 
of observed D/H ratios between different locations may be 
a too naive approach, to identify the origin of our water, 
as obviously a number of processes are at play that may 
scale differently over time for different conditions. For ex¬ 
ample, investigating molecules other than water, it has been 
found that deuterium fractionation in low-mass protostel- 
lar envelopes and molecular outflows for H 2 CO and CH 3 OH 
are higher than those for H 2 0 and NH 3 as summarized in 
Fig. 11 by |CaseH i & Ceccarelli] ( |2012[ ). This difference is ar¬ 
gued to reflect the temporal sequence in which the species 
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a b c 

Figure 1 . (a) Water molecule surrounded by four hydrogen bonded molecules as in an ideal crystal structure, (b) water molecule 
surrounded by four hydrogen bonded water molecules, where one of those contains an extra proton: ionic defect, and (c) water molecule 
surrounded by four water molecules, with one of those rotated causing a missing hydrogen bond: ‘L-defect’. 


are formed, as with decreasing gas-phase CO (increasing 
frozen out solid CO) abundances the gas-phase atomic D/H 
ratio is found to increase. Final values may also be different 
in different environments, e.g., low- vs. high-mass protostel- 
lar envelopes. In order to interpret molecular D/H ratios 
correctly, it is important to investigate all impacting pro¬ 
cesses. In the past, the focus has been on isotope dependent 
gas-phase chemical processes. So far, the two main routes 
are first gas-phase isotope dependent reactions with H 2 D + 
transferring its deuterium atom and second the enhanced 
gas-phase D/H ratio after CO freeze-out that impacts on 
ice chemistry by the higher probability of deuterium depo¬ 
sition ( Tielens|1983 Caselli fe Ceccarelli|2012 i. In the solid 
state, on the grains, deuterium fractionation is further deter¬ 
mined by the chemical surface reaction networks. These in¬ 
volve both hydrogenation and deuteration pathways and the 
crosslinks between them. This paper focuses in detail on one 
particular process only: thermally induced proton/deuteron 
exchange reactions in the ices. 

In general, the grain temperature during the formation 
of the icy mantles in the dense and dark regions of the in¬ 


voir of deuterium is HDO (Rodgers & Charnley 2002). If 


terstellar medium is low, between 10-20 K (Bergin & Tafalla 


2007). The exact temperature determines the surface reac¬ 


tion routes taken, by, e.g., sticking coefficients and the rela¬ 
tive importance of tunneling transmission coefficients, that 
lead eventually to the formation of species like CH 3 OH and 


H 2 0 ( Hiraoka et al.||1998 Watanabe fe Kouchi||2002 Fuchs 


et al. 2009 Miyauchi et al. 2008 Ioppolo et al. 20081. At 
later evolutionary stages in the star forming sequence, sur¬ 
face temperatures are no longer homogeneously distributed 


throughout the core or disk and can range up to 100 K (van 


derThkeffidJ2000 Jj0rgensenet^aLJ2OO2| Aikawaet^aLJ20d2j 


Nomura & Millar 2005| Bergin & Tafalla 2007[ Launhardt] 

et al.|2013 1. At 100 K, ices are expected to have been fully 
evaporated, thus transferring molecules from the grains to 
the gas phase. 

Within the astrochemical community, proton exchange 
in polar (i.e., water-rich) ices has attracted special attention. 
In the context of ion formation at low and higher tempera¬ 
tures, [Grim^Fal. < |1989 ) discuss the proton hopping between 
the hydrogen bonded molecules le ading to the fo rmation 
of NH)[" and OCN~. As stated by Tielens (20131 heating 


and sublimation seem to be important in interstellar ices 
near young stars in modeling the warm gas; non-energetic 
thermal reactions may therefore be relevant as well. Pro¬ 
ton exchange can be important in particular since the main 
component of ices is H a O and, moreover, the main reser¬ 


indeed such a scrambling of protons and deuterons occurs 
efficiently, this implies that not the low-temperature reac¬ 
tion routes are decisive for final HDO/H a O ratio, but rather 
the total number of H and D atoms incorporated in hydro¬ 


gen bonded molecules. Ratajczak et al. (2009) showed that 


the hydroxyl group in a methanol molecule exchanges its 
proton or deuteron with surrounding water molecules and 
proposes this as one of the possible explanations for the de¬ 
viation between modeled and observed CH 2 DOH/CH 3 OD 
ratios (Ra tajczak et al.||20 11). 

Interstellar or circumstellar ices and cometary ices 


are thought to have a common chemical origin (Bockelee- 


Morvan et al.J2000l, but the temperature processing of ices 


in comets is of an intrinsically different nature. Ices are typ¬ 
ically much thicker, i.e., orders of magnitude compared to 
the layer thicknesses of several tens of monolayers typical for 
interstellar ices. Moreover, each passage through the perihe¬ 
lion can have a large impact on the temperature depending 
on the distance to the central object. In the cometary com¬ 
munity, however, H/D exchange is often modeled only in the 


coma gas-phase chemistry via several reactions (Rodgers & 
Charnley||2002 1 following sublimation of ices. 


The mechanisms underlying proton/deuteron scram¬ 
bling have already been studied extensively in the 1980’s 
in the physical chemical community, in particular by Devlin 
and co-workers (Thornton et al.|1981 Bertie fe Devlin|1983 


|Collier et al.|1984||Wooldridge fe Devlin|1988| ) who studied 

H/D exchange in water ices as well as water-ammonia mix¬ 
tures with different isotopic compositions. For the case of 
water ice, exchange has been found to occur both in the 
bulk (Collier et al.|1984) and on the surface (Uras-Aytemiz 


et al. 2001 Park et al.|2004 Moon et al.|2010 |. In the case 

of isolated D 2 0 molecules in a H a O environment, the pre¬ 
vailing mechanism is often referred to as the so-called ‘hop- 
and-turn’ mechanism. In Fig. [T| three scenarios of the local 
ice structure are depicted: (a) a perfect crystalline structure, 
(b) an ionic defect, i.e., the presence of H 3 0 + , and (c) an 
L-defect, i.e., the lack of an expected hydrogen bond. Un¬ 
der influence of an ionic defect, the heavy water is converted 
into two coupled HDO entities that share a hydrogen bond. 
These can be converted into two nearest-neighbour HDO 
molecules by passage of an L defect. Finally, two isolated 
molecules are created after an additional proton transfer. 
Such a detailed mechanism, with corresponding distinct in¬ 
frared (IR) spectra (in the O-D stretching region) cannot be 
resolved for current astronomical IR spectra. For astrochem- 
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ical purposes, the overall reaction can therefore be summa¬ 
rized as 


Table 1 . List of (calibration) experiments and corresponding pa¬ 
rameters; temperature (T), beamline angle to the plane of the 
surface, and molecular flux (/). 


H 2 0 + D 2 0 % 2 HDO 

(Rl) 

2 HDO % D 2 0 + H 2 0. 

(R2) 


Collier et al. (19841 and Wooldridge & Devlin (19881 


studied cubic ices with D 2 0 concentrations of several per¬ 
cent in a high vacuum setup. The activation energy for pro¬ 
ton exchange was found to be ~5000 K, which means that 
the process is not likely to be relevant on interstellar time 
scales at grain temperatures up to the sublimation temper¬ 
ature of ice at 90-100 K. Since the proposed hop-and-turn 
mechanism depends on the existence of point defects present 
in the ice and on the defect concentrations, the study was 
extended to amorphous water samples by Fisher & Devlin| 
(19951. Indeed, they found a much lower activation energy 
of the turn step involved, ~3000 K. Studies of exchange on 
ice surfaces - also subject to larger structural defects - in¬ 


dicated an enhanced L-defect activity (Uras-Aytemiz et al. 
2001). These studies were performed, however, using ices 
doped with HC1 in order to have a higher sensitivity. Protons 
are generated via the reaction HC1 + H a O —> H s O + + CD, 
the exothermicity of which could influence the local pro¬ 
cesses that occur (Kim et al. 2009). Doping creates addi¬ 
tional (shallow) proton traps in the ice in the vicinity of the 
counterion ( |Uras-Aytemiz e lTal. 2001), and proton films can 
remain inactive up to around 125 K ( Lee et al.|2 007), both 
of which can influence the reaction rates detected. [Moon 
let al. |2010) showed that doping lowers the activation en¬ 
ergy found for H/D exchange on surfaces by almost a factor 
2 . 


A similar situation applies to space. In the ISM ices 
do not have a perfectly ordered structure, partially because 
they are (i) not condensed at high enough temperatures, (ii) 
comprised of many species that are formed on the surface, 
(iii) composed of more molecules than only water, and (iv) 
subject to structural changes upon energetic processing of 
the ice. These effects will contribute to a high defect con¬ 
centration in interstellar ices. 

In this paper, we aim to bring together the previous 
results of these scientific communities - physical chemistry, 
cometary chemistry and (laboratory) astrochemistry - while 
presenting several new experimental results of which the rel¬ 
evance to astronomy is discussed. We discuss amorphous 
D 2 0:H 2 0 ices with mixing ratios around 1:1. Although such 
ratios are not astrochemically relevant, they do allow for a 
higher sensitivity of the proton exchange process at lower 
temperatures, without the neccessity to use doping. We 
mimic the high expected amount of defects in the ISM by 
intentionally growing amorphous ice structures at temper¬ 
atures of 15 K and probe whether this higher defect con¬ 
centration in the ice structure, see Fig. [I leads to an in¬ 
crease of the reaction rate for reaction IRH The results are 
extended to a more general concept taking into account ‘all’ 
hydrogen bonded molecules. Subsequently, the astrochemi- 
cal implications are discussed focussing on the relevance of 
proton/deuteron swapping at long time scales, but at tem¬ 
peratures below full ice desorption. 


Experiments 
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Note. The notation a D) implies a X 10^ . 

“Instead of a full co-deposition, a layered experiment was 
performed, first co-depositing NO + H to prepare a NH 2 OH 
layer, which was then covered by D a O. ^Spectra were acquired 
at all temperatures relevant for the experiments, i.e. 90-140 K. 
“Due to the mixture preparation of HDO, the statistical ratio 
between the constituents H 2 0:HD0:D 2 0 is 1:2:1. d The effective 
deposition rate is somewhat lower as a result of HDO contami¬ 
nation in the D a O sample, this is estimated to be . 5% from the 

OH stretching mode. 


2 METHODS 

All experiments studied here consist of two sequential steps: 
(a) the simultaneous deposition of H a O (Milli-Q) and D 2 0 
(Sigma-Aldrich 99.96%) on the substrate at 15 K and subse¬ 
quently (b) an isothermal experiment at a given temperature 
(Tiao ^ 90 K) during which the level of proton exchange 
is probed by means of reflection absorption infrared spec¬ 
troscopy (RAIRS). In one experiment (Section |3.2| ) also a 
quadrupole mass spectrometer (QMS) is used. Below, the 
experimental procedure is described first, followed by an 
explanation of the calibration experiments performed, the 
analysis of the spectra, and the analysis of the temporal 
evolution of the surface abundances of H 2 0, HDO and D 2 0. 


2.1 Experimental 

Experiments are performed using the SURFRESIDE 2 setup, 
which has been constructed to systematically investigate 
solid-state reactions leading to the formation of molecules 
of astrophysical interest at cryogenic temperatures. The 
setup has been extensively described in Ioppolo et al. (2013) 
and therefore only a brief description of the procedure is 
given here. All performed experiments and the correspond¬ 
ing molecular deposition rates are listed in Table [l] The 
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Wavelength (|lm) 
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6 7 


10 




Figure 2. Bending region of the calibration spectra used for fit¬ 
ting the experiments, all spectra are recorded at 130 K. See text 
concerning the 1625 cm -1 band in the D 2 0 spectrum, the con¬ 
tamination in the HDO band, and the baseline artifact around 
1200 cm” 1 . 


Figure 3. Example of a recorded spectrum in the water bending 
region (130 K) and its best fit, decomposed into the separate 
components (H 2 0, HDO, and D 2 0 - see Fig. [2|, recorded after 
210 minutes. Note: this concerns experiment 5 from Table [T] Pure 
spectra themselves compensate for the baseline artifact at 1200 


results of experiments 1-10 are discussed in Section [3.1| Ex¬ 
periment 11 follows a procedure similar to that for the water 
experiments described below and is discussed in Section [3.2| 
The method used for experiment 12 deviates and is discussed 
in Section 13721 

SURFRESIDE 2 * 4 * * * 8 consists of three UHV chambers: one 
main chamber for ice growth, ice processing, and ice diagnos¬ 
tics, and two chambers comprising atom beam lines. All have 
a room-temperature base-pressure between 10 -9 — 10“ 0 
mbar. A rotatable gold-coated copper substrate in the center 
of the main chamber is cooled to the desired temperature us¬ 
ing a He closed-cycle cryostat with an absolute temperature 
accuracy of ^ 2 K. Both water and heavy water are pre¬ 
pared in a separate pre-pumped 10 -5 mbar) dosing line 
and after several freeze-pump-thaw cycles a co-deposition of 
room-temperature vapor H 2 0 and D 2 0 is performed. One 
species is deposited through a metal deposition line under 
an angle of 90° and the other through one of the separate 
UHV beam lines at an angle of 135° with respect to the 
plane of the surface, see Table [Tj This UHV beam line can 
be operated independently and can be separated from the 
main chamber by a metal shutter. Deposition takes place at 
a surface temperature of 15 K to ensure a large amount of 
amorphicity of the ice (see discussion above). 

After approximately 60 minutes of co-deposition at low 
temperature, yielding 45-65 ML, the substrate is heated up 
to the desired temperature, Tl so , between 90 and 140 K. 
Note that in the laboratory ices sublimate not at 100 K, 
but rather between 145-165 K, hence the higher tempera¬ 
tures employed here. The warm-up phase clearly affects the 
porosity of the ice ( |Bossa et al.||2012| and therefore also 
the ice structure. However, using H 2 0:D 2 0 mixing ratios 
around 1:1, each molecule initially has at least one neigh¬ 
bor of its isotopic counterpart, which renders the degree of 
porosity less important. Note also that the extra collapse 
of the pores at 120 K with respect to that at 90 K is only 
3% ( |Bossa eta l. 2012). Thus, assuming locally distorted hy¬ 
drogen bonded structures (Karsscmcijer et al.|2014 Karsse- 


meijer|2014 1, the exact structure of the ice and/or diffusion 

mechanisms ( |Jung et al.|2 004, Oxley et al. 2 006[ ) do not play 
a role here. A RAIR difference spectrum with respect to the 
background is acquired every 5 or 10 minutes up to the final 
time of the experiments, Use in Table [T] The background 
spectrum is acquired from an empty surface prior to the 
co-deposition at low temperature. RAIR spectra comprise a 
spectral range between 4000 and 700 cm' 1 with a spectral 
resolution of 1 cm' 1 and are averaged over 512 scans. Ex¬ 
change occurs faster at higher temperatures, therefore those 
experiments have a shorter duration. Pure component cali¬ 
bration spectra are used for the analysis of the experiments, 
see Section \2 .21 

The deposition rates mentioned in Table [I] are calcu¬ 
lated using the following relation 


cx 2 o • Rx 2 o ■ (v) 

4 • fc B • T 


(1) 


where cx 2 o is the calibration factor for the pressure gauge 
for the three isotopologues of water (X 2 0 = H a O, HDO, or 
D a O), v is the thermal velocity of the vapor molecules at 
300 K, fcn is the Boltzmann constant, and T corresponds to 
the (room) temperature. The calibration factor for water is 

1/0.9. Since there is no significant difference found between 
the absolute partial cross sections for electron-impact ion¬ 
ization of H 2 0 and D 2 0 (influencing the pressure reading of 
the gauge) ( Straub et al.|[l998[ ftikawa fe Mason||2005 1, we 
assume that this is also the case for HDO and all calibration 
factors for the pressure are taken to be equal. 

Note that the initial spectra of experiments 1-5, 7 and 

8 at 15 K after deposition are compared with each other to 
confirm the reproducibility: each low-temperature spectrum 
is analyzed with the fitting procedure described below and 
the initial surface abundances of water and deuterated water 
are found to be reproducible with a standard deviation of ~ 
5%. 
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Figure 4. Modeled H 2 0, D 2 0, and HDO concentrations for experiments 1 (a and d), 4 (b and e), and 8 (c and f) (Table |T]. The top 
row, panels (a)-(c), are zoom-ins of the HDO experimental and modeled abundance for the first 33% of the experimental duration. The 
bottom row, panels (d)-(f), depicts the evolution of all three components for the entire duration. 


2.2 Spectral Fitting 

Each spectrum is the sum of the pure H a O, HDO, and D 2 0 
components along with their intermolecular interactions. 
The intermolecular interactions cannot be captured in pure 
calibration spectra and they are stronger in the IR stretch¬ 
ing than bending region. Futhermore the bending region is 
less sensitive to crystallization effects. Therefore, our region 
of interest lies in the 2000-1000 cm -1 range, where the H 2 0, 
HDO and D a O molecules vibrate in their respective bending 
modes: 1660, 1490, 1250 cm -1 . Upon proton/deuteron ex¬ 
change the H 2 0 and D 2 0 intensities decrease, whereas the 
HDO intensity increases. 

The integrated area of each mode can be converted 
into a number of molecules, assuming that band strengths 
are available, but in reflection mode these are typically ill- 
constrained. To be able to quantify the dynamics at work, 
a spectral fitting procedure is needed to separate the three 
components. For this purpose, also three calibration exper¬ 
iments are performed. 

A calibration experiment, see final rows of Table[l] con¬ 
sists of the deposition of a ‘pure’ X 2 O component for 30 
minutes at 15 K and spectrum acquisition thereafter at all 
temperatures relevant for the experiments performed, i.e., 
90, 100, 120, 125, 130, 135, and 140 K. At each temper¬ 
ature, the spectrum looks slightly different and therefore 
pure spectra for each temperature are necessary for the fit¬ 
ting procedure. Furthermore, in the spectral region between 
~1050 and ~1250 cm -1 a known temperature-dependent ar¬ 
tifact is present, which makes it neccessary to record blank 
spectra at all relevant temperatures as well. It is most likely 
due to a temperature dependence in the absorption of our 
background sample. These spectra are also included in the 
fits (calibration experiment 4), but do not fully correct for 
the effect, as explained below. 

Using the pure components, we avoid having to correct 
for the difference in band strengths: each pure spectrum con¬ 


sists of a known amount of deposited H 2 0, HDO, and D 2 0. 
A second band around 1625 cm” 1 in the D a O spectrum has 
its origin in the combination mode lending+Vibration I Bertie 


& Whallcy 19641. There is also a small OH stretching signal 


present in the spectrum , but it is ^ 5% with respect to the 
OD stretching signal. HDO was prepared by mixing equal 
amounts of H 2 0 and D a O vapors into a pre-pumped 10” 5 
mbar) dosing line. We expect the hydrogen and deuterium 
atoms to be distributed statistically, resulting in a statistical 
1:2:1 mixture. That means that the effective HDO deposi¬ 
tion rate is only half of the total deposition rate for this 
particular experiment. The spectra shown in Fig. [2] are used 
for the fitting procedure. The HDO spectrum cannot be ob¬ 
tained pure, due to the way it is produced in the dosing 
line. In the fitting procedure we used both the original spec¬ 
trum as well as a spectrum corrected for the H z O and D 2 0 
contamination in the HDO contribution, see Section [2.4| 
Every separate experiment is composed of a time- 
resolved series of recorded spectra, each providing a snap¬ 
shot of the ongoing proton exchange process. A non-negative 
least squares solver is used to fit the pure components to the 
spectrum only in the bending range of 2000-1000 cm” 1 . Fig¬ 
ure [3] depicts an example of such a fit for the final spectrum 
(recorded after 210 min) recorded during experiment 5 from 
Table [l] In this case, the baseline does not contribute to 
the fit since the sum of the pure components themselves 
already compensate for the artifact. Usually, however, the 
temperature-dependent baseline is needed to take care of 
this artifact. 


2.3 Reaction dynamics 

The total contribution of each pure component to the fit 
is subsequently plotted versus time in order to resolve the 
dynamics. We are not able to observe the isolated HDO 
entities in a water matrix that have been proposed to be 
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produced through a hop-and-turn mechanism (Collier et al. 


19841, because both reaction partners are present in the ice 


with mixing ratios between 1:2 to 2:1. Because of this, our 
system is intrinsically easier to model, since we observe only 
the proton or deuteron hopping, leading to a reaction sys¬ 
tem of the forward and backward reactions, |R1| with rate kf 
and |R2| with rate k T . They can be evaluated in terms of a 
mean field approximation. The temporal evolutions of the 
concentrations (or abundances) can be described with a set 
of ordinary differential (rate) equations 


= —fc f [H 2 0][D 2 0] + k T [HDO] 2 (2) 

d[ ^" 0] = —fc f [H 2 0][D 2 0] + fc r [HDO] 2 (3) 

_ 2 fc f [H 3 0][D 2 0] - 2 fc r [HDO] 2 . (4) 

Where = 0. We assume thermal 

desorption to be negligable at these temperatures. Using 
an optimization algorithm in conjuction with a least-square 
ODE solver, k{ and k T can be extracted for each experiment. 
[H 2 0]o, [HDO]o, and [D 2 0]o are taken equal to the values 
found for each experiment. Section [2.4| discusses several ex¬ 
amples of both the temporal evolution of the pure compo¬ 
nents as well as the modeled result via the rate equations. 
Note that the start of each experiment is almost exclusively 
sensitive to kf, since [HDO]o is always low. We therefore 
chose to use the first 33% in time, while still in the linear 
regime, of the experiment to determine kf. A discussion on 
this topic is given in the section below. 


2.4 Optimization procedure 

Figure [4] depicts three examples of our modeled results, for 
temperatures of 90, 125, and 140 K (experiments 1, 4, and 8 
in Table [I]). The top row, panels (a)-(c), shows the HDO fit 
to the experiments for the first 33% of the experiment, and 
from the (near-) linearity of the experimental slopes it can 
be deduced that it is indeed sensitive mainly to kf. In panels 
(d)-(f) the time evolution of the experimental and modeled 
abundances (Eqns. 2)4 1 for all three species is plotted for 
the entire experimental time. 

Analyzing the abundances at 90 K we find an appar¬ 
ent discrepancy for the temporal evolution of D a O. While 


time, since there is no proton exchange at this low temper¬ 
ature, the D 2 0 abundance does not. Inspection of the IR 
spectra shows that this is likely caused by the spectroscopic 
artifact between ~1050 and ~1250 cnh 1 mentioned above. 
This artifact actually partially overlaps with the D 2 0 sig¬ 
nal. We therefore decided not to include the heavy water 
abundance explicitly in the determination of the optimal re¬ 
action rates, but to check in retrospect whether or not the 
modeled abundance is in approximate agreement with the 
experiments. 

Furthermore, the graphs at high temperature (140 K) 


somewhat from the experiment. To further constrain this, we 
have studied two limiting cases of the dependence of the H 2 0 
and D 2 0 abundances on the HDO calibration spectrum: (a) 
assuming a ‘pure’ HDO calibration spectrum and (b) assum¬ 
ing a statistical 1:2:1 distribution between H 2 0:HD0:D 2 0 


for the deposited HDO. Unfortunately, the HDO ice which 
we deposited to obtain our HDO calibration spectrum did 
not only contain HDO, but also H 2 0 and D a O. If we do 
not account for this in our calibration spectrum (case (a)), 
the water and heavy water contribution are underestimated 
during the fitting procedure. Accounting for the presence 
of both contaminants in the HDO calibration spectrum by 
assuming a statistical distribution (case (b)) will on the con¬ 
trary lead to an overestimation of the H a O and D 2 0 abun¬ 
dances. Therefore, we decided to optimize the model and 
obtain kf and k T by using only the HDO abundance. Sub¬ 
sequently, we verified that the modeled water and heavy 
water abundances fall within the range determined by the 
two limiting cases for the obtained values of kf and k T . 

Thus, the optimization sequence used for modeling the 
HDO abundance is as follows: 

(i) kf is optimized by modeling the first 33% in time of 
the experiment with k r = 0 

(ii) k r is optimized by modeling the final 33% in time of 
the experiment with kf fixed to the value found in (i) 

(iii) kf is optimized again by modeling the full experiment 
with k T fixed to the value found in (ii). 

The error function used for this optimization is 

Y' 1 ([HDO] model (t) — [HDO] exp (f )) 2 , , 

HDO 4- At <[HDO] mod ei(t) + [HDO]ex P (t)> ' 1 j 

‘■l 

The reaction rates leading to the lowest error are selected 
to be used in next iterations or are stored as the optimum 
rates. 

Finally, the activation energy of |Rl| is determined by an 
Arrhenius fit of the different reaction rates versus temper¬ 
ature. The error associated with this fit is determined with 
the stats subpackage of scipy in python and corresponds to 
the standard error of the slope. 


3 RESULTS AND DISCUSSION 

3.1 Activation energy of proton exchange in 


The reaction rates obtained from the reaction dynamical 
fitting by a simple rate equation model are converted into 
an activation energy with the use of an Arrhenius expression: 


k = v exp 


—E a 

k B T 


( 6 ) 


In this expression the prefactor, v, represents a trial fre¬ 
quency related to the vibrations of a species in a local po- 
tenti al well and correspond s to values between 10 12 — 10 13 
s _1 (Hasegawa et al. 19921. Figure [H] depicts ln(fcf) versus 


1/T. If experiments 2-10 are included in the determination 
of the activation energy, a value of 3840±125 K is found. The 
error given here is the error of the slope. The exact value 
found for the activation energy changes when the amount 
of assumed contamination in the HDO spectrum is varied, 
i. e., going from case (a) to (b). The resulting values, how¬ 
ever, are not significantly different as assessed by a two-sided 
student’s t-test (p < 0.05). 

As stated before, [Wooldridge & Devlin ( 1988| ) found a 
substantially higher activation energy of ~5000 K associated 
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Figure 5. Arrhenius plot for fcf. 


with reaction |R1| They proposed a mechanism of proton 
transfer that occurs via ionic defects in the ice. These defects 
will be more prominent in an amorphous ice and indeed the 
rate increases drastically; e - 384 °/ T / e - 500 °/ T > 4000 for T ^ 
140 K. The specific implications of this for astrochemical 
environments are discussed in Section [5] 

Experiments below 130 K are not sensitive to fc r , be¬ 
cause the situation in which enough HDO molecules are next 
to each other is never reached within the experimental dura¬ 
tion. For our experiments at 130-140 K it is not possible to 
find any correlation between fc r and the temperature, again 
because for the largest part of the experiment there are not 
enough HDO molecules close to each other. However, |Col-| 


lier et al. (1984) suggested that fc r < fcf in the case of an 


isolated D a O molecule in an H 2 0 surrounding, as a result of 
the local structure. This can also be explained from a more 
theoretical free energy point-of-view. The total Gibbs free 
energy of the reaction, G, is determined by the enthalpies, 
H, and entropies, S, of the molecules: 


AG(2 HD0-H 2 0-D 2 0) = Afff - T ■ ASf . (7) 

Using the gas-phase data reviewed in Chase| ( 1998| this leads 
to the following expression 

AG = 290 - T- 11.86 J mob 1 (8) 

= 34.88 - T- 1.43 K. (9) 

Although the thermochemistry in the solid state is affected 
by the local structure, the gas-phase expression given above 
does show that the driving force is the entropy change and 
that the reaction becomes clearly exogonic at higher temper¬ 
atures only, hence we also expect that fc r < fcf. If this relation 
is incorporated into the optimization procedure, the activa¬ 
tion energy found does not differ significantly from 3840 K. 

Finally, note that here no conclusions can be drawn 
concerning the prefactor v, because accurate values for the 
band strengths are lacking. 


3.2 Proton exchange in other hydrogen bonded 
molecules 

As mentioned in the Introduction, proton exchange in 
ices has been studied previously for several (mixtures of) 
hydrogen-bonded molecules. Here, we mention these studies 


specifically and extend the studies of water to a more general 
ice theme: proton exchange in hydrogen bonded molecules. 
As long as so-called proton wires exist between the hydrogen 


bonded molecules, exchange is expected to be rapid (Bertie 


& Devlin 1983). A proton wire can be thought of as a chain of 


hydrogen bonded molecules that has the flexibility to trans¬ 
fer protons from one to the next molecule. In such a way, any 
ionic defect in the ice can be easily passed on to the neigh¬ 
boring molecule and hence protons and deuterons can move 
efficiently. The following ice mixtures are briefly discussed: 
NH 3 :D 2 0, CD 3 0D:H 2 0, CD 3 ND 2 :H 2 0, and NH 2 0H:D 2 0. 

Exchange between isolated NH 3 in an amorphous D 2 0 
ice was found to be faster than in pure water ice and much 


faster than for crystalline ammonia (Thornton et al. 1981 
Bertie fe Devlin|1983 1. We have performed a similar experi¬ 


ment, co-depositing NH 3 and D z O (2:3) at low temperature 
and heating to 130 K, experiment 11 in Table[l] Indeed rapid 
exchange of protons can be confirmed. A steady-state HDO 
concentration is reached at least three times faster than in 
the case of water (experiment 5). This shows once more that 
amorphicity can enhance the exchange rates. 

Concerning the case of isolated amorphous d4-methanol 


mixed with water both Ratajczak et al. (2009) and Souda 


et al. (2003) confirmed thermal exchange to be efficient 


at laboratory time scales starting from 120 K. Note that 
this exchange concerns specifically the hydroxyl group of 
methanol. The C-D (or C-H) bonds do not participate in any 
hydrogen bonding network and are therefore not exchanged 
via this mechanism. They can, however, be exchanged via 
deuterium bombardment of CH 3 OH as discussed by |Na-| 
|gaoka et al. (20051. 


Ratajczak (2012) also studied the exchange of isolated 


methylamine (CD 3 ND 2 ) with water and found similar re¬ 
sults as for methanol: exchange takes place only on the 
amine entity of the molecule and is detectable from 110 K 
onwards. 

Here, in light of recent studies highlighting the hydro¬ 
genation of NO molecules and the subsequent formation 
of NH 2 OH (Congiu et al. |2012| Fedoseev et al. 20121, we 
also performed an experiment probing the exchange between 
D 2 0 and NH 2 OH (Experiment 12, Table [l]). In this case 
both hydrogen bearing groups are able to form hydrogen 
bonds. A layer of NH 2 OH is grown following a similar pro¬ 
cedure as described by Fedoseev et al. (2012) and capped 
with a layer of D 2 0 at 15 K. Increasing the sample temper¬ 
ature to 130 or 140 K could lead to exchange of protons and 
deuterons on the interface of the two molecular layers. How¬ 
ever, proton exchange is not clearly visible by RAIR spec¬ 
troscopy due to overlapping band frequencies between the 
NH(D) and OD(H) groups. Therefore, a temperature pro¬ 
grammed desorption experiment is performed, during which 
the molecules gain more mobility increasing the probability 
that NH 2 OH and D 2 0 meet, albeit for a short time until 
full desorption takes place. This experiment resulted in the 
detection of mass to charge ratios of 19, 34, 35, and 36 in the 
QMS, see Fig[6] In this Figure, the intensity does not reflect 
the in situ infrared intensity, but the intensity of gas-phase 
species after evaporation from the surface during a linear 
heating of the substrate. The temperature indicates the sur¬ 
face temperature at which this evaporation occurs. These 
masses correspond to the molecules HDO, singly and dou¬ 
bly deuterated NH 2 OH, and ND 2 OD. This confirms proton 
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Temperature (K) Temperature (K) 

Figure 6. QMS traces for TPD experiments of mass to charge ra¬ 
tios of 18, 19, 20, 33, 34, 35, and 36 probing exchange in a layered 
NH 2 0H:D 2 0 ice. Note that the x-axis changed to temperature. 


exchange in the amine and hydroxyl groups, both of which 
have hydrogen bonds. 


4 EXPERIMENTAL CONCLUSIONS 


The thermal process of proton/deuteron exchange between 
water and its isotopologues has been studied in mixed ices. 
Using RAIR spectra the time evolution of the characteristic 
vibrational bending modes of H 2 0, HDO, and D 2 0 has been 
monitered as a function of temperature. In particular, H 2 0 
and D 2 0 have been grown at 15 K and heated up to a desired 
temperature where the occurence of HDO has been probed 
indicating that the reaction H 2 0 + D 2 0 —> 2 HDO takes 
place. Temperatures between 90 and 140 K have been con¬ 
sidered and a reaction rate for each temperature has been 
determined with the use of a simple rate equation model. 
The linearity of the resulting Arrhenius plot in Fig. [5]shows 
that it is reasonable to express the rate through thermal ac¬ 
tivation and an activation energy of 3840 K has been found. 


Ratajczak (2012) suggests that proton exchange co¬ 


incides with crystallization. However, exchange has been 
found to take place prior and post crystallization, both in 
the work presented here as well as by |Wooldridge fe De-| 


vlin (19881; Fisher & Devlin (19951; Galvez et al. (20111. 


In the ISM, (F)UV radiation would render the ices largely 
compact amorphous by destruction of local molecules and 
recombination of its fragments. Exothermic energy release 
upon molecule formation is likely to have the same effect 
( Palumbo||2006 Oba et al.|2009 1. 


In this study, we mimic the lack of long-range order 
in interstellar ices by intentionally growing amorphous ice 
structures. The aim is to investigate whether a higher con¬ 
centration of structural defects (see Fig. [T|) in amorphous 
ices with respect to crystalline structures leads to exchange 
rates high enough to become relevant on interstellar time 
scales. We indeed find a lower overall (or rate-determining) 
activation energy associated with the exchange in non-doped 
amorphous ices with non-energetic detection methods com¬ 
pared to crystalline ices: 3840 versus 5000 K. This corre- 


sponds to the quantitative findings of 

Fisher & Devlin (1995|) 

for doped ices, those of Moon et al. 

2010) for pure ice sur- 


Table 2. Typical thermal exchange timescales at various grain 
temperatures and activation energies. 


T(K) 

E a = 5000 K U) 

E a = 3840 K ( 2 ) 

70 

3.3 x 10 11 yr 

2.1 x 10 4 yr 

80 

4.4 x 10 7 yr 

2.2 x 10 1 yr 

90 

4.2 x 10 4 yr 

1.1 x 10 _1 yr 


(l) 


Wooldridge & Devlin (1988 i G) This work 


faces and to the qualitative results mentioned by |Galvez| 


et al. (2011). The latter paper also discusses H/D exchange 


in water ices, but from the point of view to investigate 
whether HDO detections in ice are possible. 

Assuming a value of 10 12 s _1 for the prefactor v in 
Eqn.[6] this corresponds to typical timescales as summarized 
in Table [2] The large differences covering orders of magni¬ 
tude of the typical timescales of proton exchange clearly 
show the importance of incorporating the correct activation 
energy. 


5 ASTROCHEMICAL IMPLICATIONS 

Interstellar ices are likely to be defect-rich by nature because 
of the many ice components and various types of (energetic) 
processing. Furthermore, since ices consist mainly of water, 
the hop-and-turn mechanism is likely to dominate proton 
exchange, while, simultaneously, the existence of the proton 
wires mentioned above is key to efficient transfer. 


5.1 Protostellar and protoplanetary environments 


The typical timescales derived with the experimentally 
found activation energy can be compared to interstellar 
timescales, i.e., those mentioned in Tableland those ob¬ 
tained for protostellar and protoplanetary disk environ¬ 
ments. Schoier et al. (20021 showed that the transit time for 


grains and molecules through the warm, dense region around 
the hot core IRAS 16293-2422 is of the order of several hun¬ 
dred years. Ices present at temperatures above 80 K can 
thus be influenced by scrambling. Furthermore, although the 
high dust-temperature regions in disks concern only the in¬ 
ner few mid-plane AU ( Walsh et al.|201 0), turbulent vertical 
and radial mixing can result in transportation of water from 
the midplane to disk surface ( Furuya et al.|20l3 Albertsson 
et al.|[2014 1. This is part of a cycle in which atomic oxygen 
is transported to the midplane and reforms water and/or 
other molecules. Since the grains are transported, also the 
time they pass at higher temperatures is longer, hence al¬ 
lowing proton scrambling to take place. This then applies to 
all hydrogen bonded molecules. Mixing timescales are deter¬ 
mined by the ratio between the column density of water and 
the flux in upward or radial direction ( Furuya et al.||2013 |. 
They find that typical timescales range between 10 4 and 10' 
yrs for radii between 1 and 200 AU and that the radial ac¬ 
cretion timescale is of a similar order of magnitude and are 
long enough to allow scrambling to take place. 


Comparing typical dynamical timescales in protostellar 
and protoplanetary disk environments to the values listed in 
Ta6Ze[2[ we find that an activation energy of 3840 K renders 
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the thermally activated H/D exchange relevant at static dust 
temperatures of 70 K and above. 

We would therefore expect a D 2 0 detection probability 
lower than determined statistically in high-temperature re¬ 
gions, because any D 2 0 available in the ice is likely to be 
converted into two HDO molecules given enough time. 


In fact, consistent with this conclusion, in the warmer 
regions surrounding NGC 1333 IRAS2A |Coutens et al.| 
(20141 found a lower D 2 0/HDO ratio than in the cold en¬ 
velope layer around IRAS 16293-2422 ( Coutens et al.|2013 |. 
The higher D 2 0/HD0 than HD0/H 2 0 ratios reported in 
both cases are currently under debate. Two possible causes 
are that either (i) the surface deuteration chemistry network 
is ill-constrained or (ii) that both sublimation of grain man¬ 
tles and water formation at higher gas-phase temperatures 
takes place in the inner regions of this source. Alternatively, 
if high D 2 0 abundances compared to modeled results turn 
out to be common (Coutens et al. 2013, 2014), it must be 
because the ice mantle does not encounter high tempera¬ 
tures for long enough times. This means that the timescales 
derived here can also be used in the opposite way, to deter¬ 
mine an upper limit for the time that an icy grain resides 
in an area of a certain temperature. Currently, the limited 
number of observations does not allow us to draw such a 
conclusion. 


5.2 Cometary ices 


Shifting focus towards the application of exchange in 
cometary ices, we note that to the best of our knowledge 
no large ice chemistry models of comets exist. In the gas- 
phase coma model of Rodgers & Charnley (20021 a high ice 
abundance is sublimated as a given initial condition. The ice 
itself, however, is not modeled. They found that the D/H ra¬ 
tio in gas-phase coma species is determined by the ratio in 
their parents. Thus, the ice in fact does determine the deu¬ 
terium fractionation found in the coma. Ideally, ice abun¬ 
dances should be determined by a cometary ice chemistry 
model incorporating low- and high-temperature chemistry, 
as well as sublimation effects to investigate the effect of ac¬ 
cumulative heating on ices. 


Thermal chemistry, in particular H/D transfer, can 
play a major role in the thick cometary ices, changing the 
deuterium fractionation of many species, because HDO is the 
main deuterated component of these ices. 


The recent results on the HDO/H 2 0 ratio in the coma of 
comet 67P/Churyumov-Gerasimenko and the discrepancy 
with respect to the D/H ratios in other Jupiter family 
comets shows that it is crucial to understand the effect of 
thermal processing of ice constituents (Altwegg et al.|2014). 


groups, most of which are though to have been formed on the 
surface of dust grains, are e.g. HNO, HNC, HNCO, HNCS, 
NH 3 , nh 2 oh, CH 3 OH, HCONH 2 , CH 3 NH 2 , NH 2 CH 2 CN, 
CH 3 CH 2 OH, (CH 2 OH) 2 . 

In all examples mentioned in Section \ 0. 2\ the hydrogen 
bonded network for OH and NH moieties plays a crucial role 
in H/D exchange. Moreover, functional groups with CH that 
do not participate in hydrogen bonds do not exchange their 
proton thermally. 

We emphasize that as a first approximation exchange 
of protons and deuterons between water and other hydrogen 
bonded species should be modeled using an activation energy 
of ~ 38fO K both in disk structures as well as in cometary 
ices. 

Note that as discussed above, the forward and backward re¬ 
action rates are not expected to be equal, due to a more 
favorable entropy for the reaction H 2 0 + D 2 0 —> 2 HDO 
with respect to the backward reaction. Although the acti¬ 
vation energy itself is constrained only for the case of water 
and its deuterated analogues, arguments of enthalpy and en¬ 
tropy can also aid in constraining the reaction rates of H/D 
transfer for other molecules, as suggested by our experiments 
listed in Table |T| Finally, reactions with no net effect, such 
as H 2 0 + HDO —> HDO + H 2 0, only have meaning in a 
microscopic model, where they can be seen as analogues of 
diffusion. Thus, in a rate equation based model, they can be 
omitted. 
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5.3 Proof-of-principle modeling 

Scrambling of deuterons at high temperatures can alter the 
D/H ratios and therefore should be taken into account, when 
using D/H ratios to determine the cosmic origin of the water 
on Earth. This holds not only for the main component of 
interstellar ice, water, but especially also for any species that 
has an N-H or O-H bond that can participate in hydrogen 
bonding. Molecules detected in the ISM with such functional 
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